Neighborhood models of minority opinion spreading 
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We study the effect of finite size population in Galam's model [Eur. Phys. J. B 25 (2002) 403] 
of minority opinion spreading and introduce neighborhood models that account for local spatial 
effects. For systems of different sizes A'^, the time to reach consensus is shown to scale as In A'^ in the 
original version, while the evolution is much slower in the new neighborhood models. The threshold 
value of the initial concentration of minority supporters for the defeat of the initial majority, which 
is independent of A'' in Galam's model, goes to zero with growing system size in the neighborhood 
models. This is a consequence of the existence of a critical size for the growth of a local domain of 
minority supporters. 
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I. INTRODUCTION 

There is a growing interest among theoretical physi- 
cists in complex phenomena in fields departing from the 
traditional realm of physics research. In particular, the 
application of statistical physics methods to social phe- 
nomena is discussed in several reviews 0,0^i3, 4J. One 
of the sociological problems that attracts much attention 
is the building or the lack of consensus out of some ini- 
tial condition. There are several different models that 
simulate and analyze the dynamics of such processes in 
opinion formation, cultural d y namics , etc. [3 , B, B , B , 
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Among all those models, the one introduced by Galam 
[3, i^i to describe the spreading of a minority opinion, in- 
corporates basic mechanisms of social inertia, resulting 
in democratic rejection of social reforms initially favored 
by a majority. In this model, individuals gather during 
their social life in meeting cells of different sizes where 
they discuss about a topic until a final decision, in favor 
or against, is taken by the entire group. The decision is 
based on the majority rule such that everybody in the 
meeting cell adopts the opinion of the majority. Galam 
introduced the idea of "social inertia" in the form of a 
bias corresponding to a resistance to changes or reforms, 
that is: in case of tie, one of the decisions (in the orig- 
inal version, the one against) is systematically adopted. 
We will describe in detail the model and its main con- 
clusions in the next sections. This simple model is able 
to explain why an initially minority opinion can become 
a majority in the long run. An interesting example was 
its application to the spread of rumors concerning some 
September 11-th opinions in France 8]. One of the ma- 
jor conclusions of the mean-field-like analysis in Ref.Q, is 
the existence of a threshold value Pc < 1/2 for the initial 
concentration of individuals with the minority opinion 
(against the social reform). For p > pc every individual 
eventually adopts the opinion of the initial minority, so 
that the social reform is rejected and the status quo is 
maintained. A related message of this result is that a ru- 
mor spreads, although initially supported by a minority, 
if the society has some bias towards accepting it. 



Galam traces back his results to dynamical effects pro- 
duced by the existence of asymmetric unstable points 
previously considered by Granovetter and Schelling 
j24j . These are fixed points of recursion relations describ- 
ing the dynamics of the fraction of a population adopting 
one of two possible choices. In threshold models these re- 
lations are obtained considering a mean field type of in- 
teraction in which individual thresholds to change choice 
(tolerance) are compared with the fraction of the popula- 
tion that has already adopted the new choice. Granovet- 
ter himself discusses that the stability of the fixed 
points can be changed by spatial effects, noting that the 
assumption that each individual is responsive to the be- 
havior of all the others is often inappropriate. In such 
complete connectedness of the population seems to be cir- 
cumvented by the introduction of the meeting cells. Only 
individuals in each meeting cell interact among them- 
selves. In this sense each meeting cell plays the role of 
a bounded neighborhood and it is still possible to 
obtain analytically recursion relations for the dynamics. 
However, contrary to the bounded neighborhood model 
of Schelling, individuals enter and leave these neighbor- 
hoods randomly, and the neighborhoods do not have any 
characteristic identity other than their sizes. Even if the 
meeting cells are thought of as sites where local discus- 
sions take place, Galam's model does not incorporate 
local interactions since the individuals are randomly re- 
distributed in the meeting cells at each time step of the 
dynamics. 

The alternative considered by Schelling to the bounded 
neighborhood model is a spatial proximity model in 
which everybody defines his neighborhood by reference 
to his own spatial location. The spatial arrangement 
or configuration within the neighborhood mediates the 
interactions. We propose here a different neighbor- 
hood model which shares some characteristics with the 
bounded neighborhood and spatial proximity models: 
The meeting cells are neighborhoods defined by spatial 
location, therefore introducing important local effects, 
but the interaction within the neighborhood is indepen- 
dent of the spatial configuration within the cell. Con- 
trary to the model in the individuals are here located 
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at fixed sites of a lattice. The local neighborhood or 
meeting cell in which a given individual interacts changes 
with time, reflecting neighborhoods of changing shape 
and size. Such neighborhood model could be appropri- 
ate for a relatively primitive society in which interactions 
are predominantly among neighbors, but the size of the 
neighborhood or interaction range is not fixed. 

A different version of Galam's model was introduced by 
Stauffer Q . At variance with our neighborhood models 
in which individuals are fixed in the sites of a lattice and 
the meeting cells have a maximum size, Stauffer consid- 
ers the situation in which individuals freely diffuse in a 
lattice with only a fraction of sites being occupied. This 
diffusion process leads to the formation of "natural" clus- 
ters which play the role of our meeting cells: It is within 
each one of these clusters in which the rule of major- 
ity opinion and bias towards minority in case of a tie are 
taken. Stauffer finds in his model that the time to reach a 
consensus opinion grows logarithmically with system size 
N. We also find this dependence in the original model of 
Galam, while in our neighborhood models the consensus 
time takes much larger values and is compatible with a 
power law dependence. A related model, including the 
figure of the "contrarians" (that is, people that always 
oppose to the majority position), was later introduced 
by Galam 25] and also analyzed by Stauffer "2^. 

A main consequence of introducing the spatial ef- 
fects considered in our neighborhood models is that the 
threshold found in 0] disappears with system size, i.e. 
limAT^co Pc = 0, so that in large systems the minority 
opinion always spreads and overcomes the initial major- 
ity for whatever initial proportion of the minority opin- 
ion. This is a consequence of the existence of a critical 
size for a local domain of minority supporters. Domains 
of size larger than the critical one will expand and oc- 
cupy the whole system. For large systems there is al- 
ways a finite probability to have a domain of over-critical 
size in the initial condition. While in traditional Statisti- 
cal Physics we are mostly concerned with the thermody- 
namic limit of large systems, these findings emphasize the 
important role of system size in the sociological context 
of models of interacting individual entities. 

The outline of the paper is as follows. Section 2 re- 
views the original definition of Galam's model 0,01 and 
introduces our new local neighborhood models. In Sect. 

3 we go beyond the mean field hmit of Refs.0,11] by dis- 
cussing the system size dependence of the predictions of 
the original model. Steady-state and dynamical proper- 
ties of our neighborhood models are presented in Sects. 

4 and 5. General conclusions are summarized in Sect. 6. 



II. DEFINITION OF THE MODEL 



defined by the number of individuals k that can meet 
in the cell. Let us define Uk as the probability that a 
particular person is found in a cell of size k. Obviously, 

it is J2k «fc = 1- 

The dynamics of the model is as follows: first the meet- 
ing cells are defined by giving each one a size according 
to the probability distribution {a^}, such that the sum 
of all the cell sizes equals the number of individuals N, 
but otherwise their location or shape are not specified. 
These cells are not modified during the whole dynamical 
process. The persons have an initial binary (against, 
or in favor, — ) opinion about a certain topic. The prob- 
ability that a person shares the -I- opinion at time t is 
P+{t) and an equivalent definition for P-{t) = 1 — P+{t). 
Initially one sets P+{t = 0) = p. Alternatively, P+{t) 
can be thought of as the proportion of people supporting 
opinion + at time t. 

The N individuals are then distributed randomly 
among the different cells. The basic premise of the model 
is that all the people within a cell adopt the opinion of 
the majority of the cell. Furthermore, in the case of a 
tie (which can only occur if the cell size k is an even 
number), one of the opinions, that we arbitrarily identify 
with the -|- opinion, is adopted. Once an opinion within 
the different cells has been taken, time increases by one, 
t ^ t + 1 and the individuals rearrange by distributing 
themselves again randomly among the different cells. 

The main finding of this model is that an initially mi- 
nority opinion, corresponding to p < 1/2 can win in the 
long term. This is an effect of the tie rule that selects 
the -t- opinion in case of a tie. 

It is possible to write down a recursion relation for the 
density of people that at time t have the -I- opinion as 

/Ic\ 
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the notation [x] indicates the integer part of x. This is a 
mean-field equation that neglects possible fiuctuations. 

For a wide range of distributions {afc} this map has 
three fixed points: two stable ones at P+ — 1 and P+ = 
and an unstable one, the faith point, at P+ = pc- Hence, 
the dynamics is such that 



lim P+{t) = 



if 
if 



> Pc 
< Pc 



(3) 



A. Galam's original non-local model 



B. Neighborhood models 



The model considers a population of N individuals who 
randomly gather in "meeting cells" . A meeting cell is just 



We now introduce our Neighborhood Models that in- 
corporate local spatial effects in the interacting dynamics 
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FIG. 1: (a) Regular 2D tessellation; All the cells are simul- 
taneously created, being tUx and my uniformly distributed 
between 1 and M. (b) Locally grown tessellation: A site 
is chosen, and from it a cell of size mx,my, excluding 
those already belonging to other cell. 

proposed by Galam. In these local models, individuals 
are fixed at the sites of a regular lattice and they interact 
with other individuals in their spatial neighborhood. We 
have considered several cases: 



1. One- dimensional neighborhood model: synchronous 
update 

The N individuals are distributed at the sites of a lin- 
ear lattice i = 1,2, . . . ,N . Once distributed, they never 
move again. Initially they are assigned a probability p of 
adopting the -I- opinion, and (1 — p) the — opinion. The 
dynamics starts by defining the meeting cells /c = 1, 2, . . . 
as the segments [ik,ik+i ~ 1] of length Wfc = ik+i — ik- 
The cell sizes nik are distributed according to a uniform 
distribution in the interval [1,M]. The average cell size 

is hence (mk) = ^^^2*'^'' and the average number of cells 
is N/{mk)- Once the cells are defined, the dynamical 
rules of Galam's model are applied synchronously to all 
the cells, time increases by one t ^ t + 1. In the next 



time step, new cells, uncorrelated to the previous ones 
are defined and the dynamical rules applied again. The 
process continues until there is a single common opinion 
in the whole system. 



2. Two-dimensional neighborhood model: synchronous 
update 

The two-dimensional case is very similar to the ID ver- 
sion explained above. The only difference is the way the 
meeting cells are defined in the two-dimensional lattice. 
An individual is now characterized by two indexes 
with I < i, j < L, such that the total number of indi- 
viduals is iV = L^. We have considered two different 
definitions of the cells originated in two tessellations of 
the plane: (a) the regu7ar tessellation and (b) the locally- 
grown tessellation. In the regular tessellation, we de- 
fine segments in the i and j axis independently, such 
that the sizes are in both cases uniformly distributed be- 
tween 1 and M in the same way that we did in the one- 
dimensional case. Figure^ plots a typical example. In 
the locally-grown tessellation we first choose a site of the 
lattice and then define a rectangle around it whose sides 
are both uniformly distributed between 1 and M. The 
cell is defined then as the sites in the resulting rectangle 
excluding those sites that already were part of a previ- 
ously defined cell. Figure QJ) shows a typical example. 

Once the cells are defined, the dynamical rules are ap- 
plied synchronously to all the cells and a common opinion 
is formed within each cell. Time then increases by one 
t —f t-\-l. In the next time step, new cells are defined and 
the process continues until a consensus opinion is reached 
in the whole population. 

3. Asynchronous update 

The ID and 2D models have been also considered in 
the asynchronous update version. In this case, a lat- 
tice site is randomly chosen and a cell defined around it 
as a segment (ID) of size to or a rectangle (2D) of size 
{mx,my). It is only within this cell that the biased ma- 
jority rule is applied. Time increases hy t — > t -\- m/N 
in ID and hy t t + {mj.my)/N in 2D. Then a new 
site is selected randomly and the process iterates until a 
consensus opinion is obtained. 



III. RESULTS FOR GALAM'S ORIGINAL 
NON-LOCAL MODEL. 

We present in this section an analysis of Galam's orig- 
inal non-local model. Our aim is to go beyond the mean- 
field approach of references 0, IE by studying the 
system size dependence of the different magnitudes of 
interest. Some of the results are based on numerical sim- 
ulations of the model. 



4 



80 



60 



§40 



20 



1 1 1 1 1 1 

consensus to 


;l ' 1 




- 




consensus to 


I.I.I 


+ 

I.I 



I I I , I , I , I , L- 

1.1 0.2 0.3 0.4 0.5 0.6 



FIG. 2: Phase diagram of the original Galam's model in the 
plane (pc, M). 



We consider N individuals that distribute themselves 
randomly in meeting cells whose size is uniformly dis- 
tributed between 1 and M. In the notation of Eq. (1), 
this means that at = 2fc/(M(Af + 1)), fc G [l,Af], as it 
follows from the fact that Uk measures the probability 
for an individual of being in any of the k sites of a cell 
of size k. Initially we assign to each of the persons any 
of the two possible opinions, such that the probability of 
having the favored opinion is p. Again, in the language 
of Eq. (P), we are setting P+(0) = p. We then apply 
the dynamical rules of Galam's model until a consensus 
opinion is formed. By iteration of this procedure, we 
measure the probability p that the consensus opinion co- 
incides with the favored one, +. This is precisely defined 
as the fraction of realizations that end up in the favored 
+ opinion. 

The analysis of Eq. predicts a first order phase 
transition in the sense that the "order parameter" p = Q 
if p < Pc and p = 1 \i p > Pc- In Fig. |21we show, in 
the parameter space {pctM), the regions where the two 
solutions, as obtained by finding numerically the non- 
trivial fixed point of the recurrence Eq. exists. Note 
that, as expected, the larger the decision cells, the closer 
the faith point to 1/2. Notice also in this figure that 
some pairs of consecutive values of M give almost the 
same value for pc- The reason is that, for odd values of 
M, the rule that applies in case of a tie is used only up 
to the value M — 1 (which is even), so odd numbers give 
similar values of pc than the precedent even number. 

Since Eq. ^ neglects possible system size fiuctuations, 
Eq. (PI this result is only valid in the limit N ^ oo. We 
plot in Fig. the raw results of our simulations for 
different system sizes. The analysis carried out in Fig. 
|3Jd shows that the asymptotic results of iV cx) are 




(P-PJN' 

FIG. 3: (a) Order parameter in Galam's non-local model. 
The white symbols correspond to the case Al — 4 {pc = 
0.2133077908...), while the black ones to M = 5 (pc = 
0.3467871056 . . .). The values of range between TV = 10^ 
and A'^ = 10^. (b) The order parameter for both values of M 
for rescaled values of p according to (p — Pc)N^^'^. Also, it is 
shown a fit with the function p{p,N) = (1 + erf(a;/1.17))/2 
with the scaling variable x = {p — pc)/N^'"^ . 

achieved by means of the following scaling relation 

p(p,7V) = /((p-Pe)iV-i/2). (4) 

Therefore, there is a region of size N^^^^ where there is 
a significant probability that the results differ from the 
infinite size limit. 

We now analyze the time T it takes to reach the consen- 
sus opinion. Strictly speaking, in the mean field approach 
the number of iterations needed for Eq. (1) to converge 
to the fixed point is infinity. In Ref. it was adopted 
the criterion that the fixed point had been reached at the 



5 



40 r 



30 



t-^ 20 



10 




0.2 



0.3 
P 



0.4 



0.5 




FIG. 4: Time that it takes to the system to reach the consen- 
sus state as a function of the initial probability p for differ- 
ent system sizes. The white symbols correspond to the case 
M = 4 while the black ones to M = 5. The values of A'' range 
between = 10^ and iV = 10^ 



FIG. 5: Time to reach the consensus state as a function of 
system size for different values of p and M = 4. Two different 
slopes can be seen for values of p greater and lower than the 
critical value pc ~ 0.21. Also the predicted slopes are plotted. 



time T such that P+{T) differs from the fixed point in 
less that 0.01. In our simulations, and as in reference [T3| . 
it is natural to define the time T as the finite number of 
steps needed to reach the consensus opinion. We plot in 
Fig. ^the time T as a function of the initial probability 
p of the favored opinion. It can be seen that the time T 
increases with increasing system size and takes its maxi- 
mum value at the faith point. A closer analysis shown in 
figure shows that, for all values of p, the time needed 
to reach the consensus increases logarithmically with N . 
It is interesting to note that the same logarithmic depen- 
dence was found in Ref. 0] for the model accounting 
for Brownian diffusion. 

A simple argument can help to understand this loga- 
rithmic behavior. One can mimic the size dependence of 
the time T by noticing that the definition used in the nu- 
merical simulations is equivalent to define T as the value 
for which P+(T) = 0{l/N) or 1 - P+{T) = 0{l/N), 
since 1 /N is the minimum possible value for p over p — Q. 
This can be now obtained by linearizing the evolution 
equation around any of the fixed points 0, 1. Defining re- 
spectively 5{t) — P+(t)— or 5{t) — P+(t) — 1, and replac- 
ing in the evolution equation Q P+{t + 1) = F[P+{t)\: 
we obtain at first order: 

6{t + l) = \5{t), (5) 

where A = F'[0] and A = f [1], respectively. In this linear 
regime, we have simply: 5{t) = A* 5(0) and according to 
the definition above. 

In 5(0) -In TV 

InA ' ^ ' 



where the correct value of A has to be used in each case. 
This is the logarithmic behavior observed in the simula- 
tions. Furthermore, the slope of the logarithmic law will 
depend on the fixed point at which the system tends. 
So, at one side of the critical point, the slope should be 
different than at the other. Figure |S1 shows the confir- 
mation of this prediction, where it is seen that for high 
values of p (above the faith point) the slope is quite dif- 
ferent than for lower values. For M = 4, the predicted 
values for A are A = 1/10 and A = 1/2 for the fixed point 
at P+ = 1 and P+ = respectively. The corresponding 
slopes, -1/lnA, are 0.434... and 1.442.... As shown 
in the figure, these values agree well with the measured 
slopes. The only discrepancy is for p ^ Pc for which the 
time needed to reach consensus must include as well the 
time needed to leave the fixed point Pc- 



IV. NEIGHBORHOOD MODELS: STEADY 
STATE PROPERTIES 

In this section we consider the steady state properties 
of our neighborhood models defined in Sect. IIIBI We will 
see that the introduction of spatial local effects leads to 
a very different behavior than for the non-local version of 
the previous section in which individuals were distributed 
randomly in fixed cells. 

In our neighborhood versions, similarly to the origi- 
nal Galam model, it turns out that a consensus opinion 
is always reached in a finite number of steps. We first 
consider the order parameter, p, defined as the proba- 
bility that the consensus opinion coincides with the fa- 
vored one. Figures and [3 show that, both in the ID 
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FIG. 6: Order parameter for the neighborhood one- 
dimensional system with synchronous update for M — 5. The 
system size ranges between A'^ = 10 (rightmost curve) and 
N = 10"* (leftmost curve) . The inset shows the validity of the 
scaling law p = p{pN") with a = 0.7. 
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FIG. 7; Same as Fig. El for neighborhood models with syn- 
chronous update and regular 2D tessellation. M — 5. The 
system size is N = L x L with L between L — 15 (rightmost 
curve) and L = 10"^ (leftmost curve). The inset shows the 
scaling law with a = 0.24. 



and the 2D cases, the order parameter p is an increas- 
ing function of the initial probability p of adopting the 
favored opinion. It is possible to define, quite arbitrar- 
ily, the transition point pc as the one for which p = 1/2. 
However, pc depends upon the system size as a power 
law Pc{N) ~ N~°', hence for increasing N the transition 
point tends to Pc = 0. In other words, the transition dis- 
appears in the thermodynamic hmit, N — > oo, and the 
favored opinion, in that limit, is always the selected one 
independently of the initial choice. More precisely, the 
data can be described by the scaling law 

p{p,N)=pipN^). (7) 

The exponent a depends on the dimension and on the 
maximum size M of the domains. In Fig. |Slwe plot the 
Af-dependence of a in several ID and 2D neighborhood 
models. Notice that a decreases for increasing M and 
it depends on the system dimension, but it is basically 
independent of the local rules defined. Alternatively, for 
fixed p we can define a critical value Nc{p) such that a 
small population N < Nc{p) tends not to propagate the 
initially minority opinion. 

V. NEIGHBORHOOD MODELS: DYNAMICAL 
EVOLUTION 

There are several differences between the evolution in 
the neighborhood and non-local models. We first ana- 
lyze the time T needed to reach the consensus. As in 
Galam's original model, the data in Figs. El and ^1 for 
the ID and 2D cases, respectively, show that for fixed 
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FIG. 8: Scaling exponent defined in equation Q as a 
function of A/"*. The different curves correspond to: ID 
(one-dimensional system, synchronous update), while the bi- 
dimensional are; 2D-1 (synchronous update, regular tessella- 
tion), 2D-2 (synchronous update, locally grown tessellation), 
2D-3 (asynchronous update). 

N, the time T reaches a maximum at the critical point 
Pc{N). Notice that the numerical values for T are much 
larger in the local models that they were in the original 
model. Furthermore, as shown in the insets of Figs. El 
and ^1 for fixed p, the time T has two different growth 
laws according to whether the population N is smaller 
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FIG. 11: Plot of the dynamical evolution of a non-local 
system updated according to the synchronous update. Each 
frame correspond to a time, increasing from left to right and 
top to bottom. The system size is iV = 256^ and M = 5. 



FIG. 9: Time to reach the consensus vs p, for synchronous 
ID local cells and for M = 5. Same symbols and systems sizes 
than in Fig. El The inset shows the time to reach consensus 
plotted against system size, for different values of p, i.e. p — 
0.2 (triangles), p = 0.1 (diamonds), p = 0.05 (circles), p = 
0.02 (squares). Two power- laws are observed (depending on 
whether TV < iVc(p) or iV > iVc(p)); for N < Nc, P ^ 1.6 and 
in the regime N > Nc{p), P ^ 0.2. 




FIG. 10: Time to reach the final consensus state vs p, for 
the 2D cells, and for M = 5. Same symbols and systems sizes 
than in Fig. The inset shows the time to reach consensus 
for fixed values of p: p = 0.3 (stars), p — 0.2 (triangles), p = 
0.1 (diamonds), p — 0.05 (circles), p = 0.02 (squares). The 
results are quite similar to those depicted in Figure |5] In this 
case, for N < Nc{p), /3 ~ 0.6 and in the regime A'^ > Nc{p) (a 
regime only clearly seen in the cases p = 0.3, 0.2), /? ^ 0.1. 



or larger than the critical size Nc{p). For N < Nc{p) the 
time T increases as a power-law T ~ of the system 
size N with [3 = 1.6(0.6) for d = 1(2). For iV > Nc{p) 
the data are compatible with a power law with smaller 
exponents, although it can not be completely excluded a 
logarithmic dependence in this case. 

In Fig. ^2 we plot several snapshots corresponding 
to the evolution according to Galam's original rules. It 
is seen that the evolution is very fast and at each time 
step the number of people favoring a particular opinion 
increases but due to the non- locality of the rules, one can 
not see any structured pattern of growth. 

When using the neighborhood rules, however, it can be 
seen that, after a very short transient in which domains 
are formed, the ulterior evolution is by modification of 
the interfaces between the two possible types of domains. 

Figure IT^ depicts the dynamical evolution of the syn- 
chronous ID version of the neighborhood model. It shows 
a linear growth of the size of domains of favored opinion. 
It is easy to predict the slope of the linear growth ac- 
cording to the following reasoning: the only evolution 
is produced if the discussion cell intercepts two domains 
with different opinions. Let us consider the position of 
an interface at time t = Q. For cells of odd size, the bias 
rule does not apply and the interface does not move in 
the average. For cells of even size, the bias rule acts only 
in one case and it is easy to show that the location of the 
interface, on the average, increases by 1/2. Therefore, 
the size of the domain of the favored opinion increase as 
t/2. Averaging for an equiprobability of having cells of 
size between 1 an M, yields a linear growth (it with 



/3 = 2 



[f ] [f + 1] 
M{M + 1) 



(8) 



where [a;] once again denotes the integer part of x. The 
validity of this result for M = 5 is shown in Fig. 1121 

The 2D snapshots of the dynamical evolution, see Fig. 
1131 show that as in the ID version, domains are formed in 
an initial transient, and after this, grow by interface dy- 
namics. To characterize the growth of the characteristic 
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FIG. 12: Plot of the dynamical evolution of a local ID system. 
Each row corresponds to a time, increasing from bottom to 
top. The system size is A'' = 5000 and M = 5. It is also shown 
the predicted slope given by Eq. JHJ The inset amplifies a 
small region for a short time evolution. 



size of the domains we calculate a characteristic radius 
Rc defined as 



J2^ ^+{iyi 



(9) 



Where S^{i) is 1 if individual i supports opinion +, 
otherwise and is the position in the square lattice of 
th individual. Figure FPil shows that this characteristic 
linear dimension grows linearly with time as in the ID 
model. 

The growth of domains in ID and 2D occurs for those 
domains whose initial characteristic size is larger than 
a critical one R*. When the initial size is larger than 
i?*, domains grow until they take over the whole system, 
while smaller domains shrink toward extinction. The 
quantitative calculation of the critical radius starts by 
placing a unique circular island of radius r of minority 
people surrounded by the majority and watching it grow 
or shrink. This initial condition is quite different of that 
of a random initial distribution. In order to avoid finite 
size effects, one has to be careful to choose a system size 
A'' such that the initial density p = Tir^ /N is lower than 
the critical one, Pc- In the inset of Fig. ^| we show 
that the probability that an initial circle of radius 

r grows, is indeed independent of system size. We define 
the critical radius as the value of r such that p{r) = 0.5. 
Figure [131 shows that the value of the critical radius in- 
creases linearly with ill, both in ID and in 2D models. 




FIG. 13: The dynamical evolution of a 2D local model up- 
dated according to the asynchronous rule is shown. In the 
simulation iV = 400 x 400, M = 5 and p = 0.2. The time 
increases from left to right and from top to bottom. The 
snapshots correspond, respectively to T = 0, 20, 52, 106, 276, 
394, 594 and 818. 




FIG. 14: Characteristic radius, Rc averaged over 100 in- 
dependent runs as a function of time for a two-dimensional 
system updated according to the synchronous update and reg- 
ular tessellation. Rc is shown to grow linearly with time for 
different values of M . The initial condition is taken as a single 
circular domain of minority supporters with a radius R> R* . 



VI. GENERAL CONCLUSIONS 

We have revisited Galam's model |3i Hi of minority 
opinion spreading and introduced related neighborhood 
models that incorporate spatial local effects in the inter- 
actions. These models share basic characteristics with 
the bounded neighborhood and spatial proximity mod- 



els of Schelling In both cases we have considered 

in detail the role of system size in the properties of the 
system. For the original nonlocal version T, , we have 
found that the transition from initial minority final dom- 
inance to initial minority disappearance is smeared out 
in a region of size N~^/'^, while the time it takes to reach 
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FIG. 15: Critical radius vs M, for the different versions of 
the neighborhood systems: (a) ID model, (b) 2D-1, 2D-2, 
2D-3 as defined in Fig. |H| The insets show the probability 
p(r) that an initial domain of radius r grows; (a) corresponds 
to M = 15 and (b) is for the 2D-1 model for M = 7, both for 
different system sizes. 



complete consensus increases as IniV, as in StaufFer's per- 
colation model [l^ . 

In our local neighborhood models, we have considered 
ID and 2D lattices with regular and locally grown tes- 
sellations, both with synchronous and asynchronous up- 
dates. All these local versions behave qualitatively in 
the same way. The most important finding is that the 
threshold value for the initial minority concentration pc 



decreases as iV~", such that the transition from initial 
minority spreading to majority dominance disappears in 
the thermodynamic limit N ^ oo. The neighborhood 
models are, in this sense, more efficient to spread an ini- 
tially minority opinion. However, while a nonlocal model 
is very fast in spreading a rumor, the corresponding re- 
laxation times to reach consensus in the neighborhood 
models are much larger since they turn out to increase 
with a power of the system size TV. 

We have also shown that the fact that liniN^ooPc = 
is due to the existence of a critical size for a spatial do- 
main of minority supporters. For large enough systems 
there is always an over-critical domain that spreads and 
occupies the whole system, with a characteristic aver- 
age radius growing linearly with time. This critical size 
domain has some analogies with the critical nucleus of 
nucleation theory 27], but it has different characteris- 
tics. In classical nucleation the existence of a critical 
radius is due to the competition between surface tension 
and different bulk energy between the two possible ho- 
mogeneous states. The concept of critical nucleus in this 
context is not meaningful for one-dimensional systems for 
which no surface tension exists. An over-critical droplet 
appears as a rare fluctuation in the bulk of a metastable 
state and it then grows deterministically. Noisy pertur- 
bations in the growth dynamics are generally a second 
order effect. In our case over-critical domains appear in 
the random initial condition. The critical size is here 
an average concept resulting from the competition of the 
bias favoring the minority opinion in case of a tie (analog 
of bulk energy difference) and a stochastic dynamics that 
might lead to the disappearance of the minority domain, 
surface tension being a second order effect. Critical size 
means here equal probability for the domain to spread or 
to collapse. 

The existence of this critical size clearly shows an 
important difference between typical statistical physics 
problems and sociophysical ones. In the former case, 
one is mostly concerned with the thermodynamic limit 
of large systems, while these findings emphasize the im- 
portant role of system size in the latter case. 
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